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Abstract. Most stars - especially young stars - are observed to be in multiple systems. Dynamical evolution is unable to 
pair stars efficiently, which leads to the conclusion that star-forming cores must usually fragment into > 2 stars. However, the 
dynamical decay of systems with > 3 or 4 stars would result in a large single-star population that is not seen in the young stellar 
population. Additionally, ejections would produce a significant population of close binaries that are not observed. This leads to 
a strong constraint on star formation theories that cores must typically produce only 2 or 3 stars. This conclusion is in sharp 
disagreement with the results of currently available numerical simulations that follow the fragmentation of molecular cores 
and typically predict the formation of 5-10 seeds per core. In addition, open cluster remnants may account for the majority of 
observed highly hierarchical higher-order multiple systems in the field. 
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1. Introduction 

In the Galactic field most solar-type stars are in multiple sys- 
tems with a multiplicity fraction, f mu \ t = (B + T + Q)/(S + B + 
T + Q) — 0.58 where S,B,T,Q are the number of single, bi- 
nary, triple and quadruple systems, respectively (Duquennoy 
& Mayor 1991, including significant corrections for incom- 
pleteness from the raw value of 0.49). Recent studies have 
found that / mu i t in the field may be somewhat higher than 
this (Quist & Lindegren 2000; Soderhjelm 2000). Duquennoy 
& Mayor (1991) find the uncorrected ratio of systems to be 
S : B : T : Q = 1.28 : 1 : 0.175 : 0.05 (see also Tokovinin & 
Smekhov 2002). 

There is evidence that / mu i t is significantly higher among 
young stars than in the field (e.g. Leinert et al. 1993; Mathieu 
1994; Duchene 1999; Bodenheimer et al. 2000; Patience et al. 
2002; Kroupa et al. 2003; Duchene et al. 2004; Haisch et al. 
2004). In Taurus essentially all stars between 0.3 and 1 M Q 
are in multiple systems with separations extending to about 
1500 au (Leinert et al. 1993; Duchene et al. 2004; Haisch et al. 
2004). Samples of T Tauri stars in other nearby star-forming 
regions also show that T Tauri stars generally have a high / mu i t 
and a bias towards larger separations than the main-sequence 
sample (Duchene 1999; Patience et al. 2002; Duchene et al. 
2004; Haisch et al. 2004). Such samples are generally biased 
towards large separations (~ 100 to 1000 - 3000 au) and may 
miss many close binaries. Generally it is found that / mu i t in the 
sampled separation ranges is approximately twice that in the 
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Galactic field and if extrapolated to the whole separation range 
implies / mult ~ 1. 

It has been shown that it is not possible to reproduce the 
observed / mu i t through the dynamical evolution of star clus- 
ters that are born with a single-star population (Kroupa 1995a). 
Dynamical interactions are able to disrupt many wide binaries, 
but are not able to pair stars efficiently or significantly change 
the properties of close binaries (e.g. Kroupa 1995a,b; Kroupa 
et al. 1999; Kroupa & Burkert 2001; Kroupa et al. 2001). This 
leads us to the conclusion that the majority of stars must form 
in multiple systems. 

All stars form in dense molecular cores (Andre et al. 
2000), and most of these star-forming cores are found within 
larger complexes that will become star clusters. It is these 
star-forming cores that must fragment into multiple objects. 
Recent simulations have shown that cores are able to fragment 
if there is a small amount of initial rotation (eg. Burkert & 
Bodenheimer 1996), with compression due to an external pres- 
sure (eg. Hennebelle et al. 2003, 2004), or if they contain high- 
(eg. Bate et al. 2002, 2003; Delgado Donate et al. 2004) or low- 
(eg. Goodwin et al. 2004a,b) levels of turbulence. As reviewed 
by Duchene et al. (2004) the models predict the cores to gener- 
ally fragment into N = 5 - 10 pre-stellar objects. But Duchene 
et al. point out that this appears to contradict the observational 
findings according to which no, or very little, higher-order sub- 
clustering within cores is evident. 

In this paper we re-address this problem by considering the 
elementary constraints on the number of stars that can form 
within a core posed by the decay of higher-order multiple sys- 
tems. Limits on the number of objects that usually form within 
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cores provide very strong constraints on theories of star forma- 
tion. 

2. Constraints on the fragmentation of cores 

As a core collapses it is initially isothermal as it is optically thin 
to its own radiation. At a critical density of p cr { t ~ 10~ 13 g crrT 3 
for a 1 M & core, the gas becomes optically thick and enters an 
adiabatic regime (eg. Larson 1969; Tohline 1982; Masunaga & 
Inutsuka 2000). This results in a minimum Jeans mass being 
reached around p cr ; t of M m j n - 10~ 2 M©. It is at ~ p cr ; t that frag- 
mentation is expected to occur producing fragments of mass 
M m ; n . Thus the length scale of the initial multiple systems is 
expected to be 

(1 A/I \ 1/3 
< 125(M core /M ) 1/3 au, (1) 
47TPcrit ) 

where M core is the mass of the pre-stellar core. This scenario is 
supported by the observed peak in the T Tauri binary distribu- 
tion at ~ 100 au (eg. Mathieu 1994; Patience et al. 2002). 

2.1. Dynamical decay 

Any (non-hierarchical) multiple system with N > 3 members 
is generally unstable to decay with a half-life of fdecay ~ 80 
crossing times, 

fde -~ 14 U UrJ yrs ' (2) 

where R is the size of the system and M stars the mass of the com- 
ponents (Anosova 1986; see also Reipurth & Clarke 2001). A 
typical decay time-scale for a system with R = 250 au (2Rf olm ) 
and M stars = 1M is then of order 55 kyr. The decay of mul- 
tiple systems has two important consequences: the preferential 
ejection of the lowest-mass member of the system (Anosova 
1986; Sterzik & Durisen 2003), and the reduction of the semi- 
major axis of the remaining binary by a factor of 2-4 (Anosova 
1986). The reduction in the semi-major axis is significantly 
greater if ejections occur whilst accretion is on-going (Umbreit 
et al. 2005). Thus, for decaying higher-order systems fdecay falls 
rapidly as each ejection reduces the size of the system signifi- 
cantly (as fdecay x R 3 ^ 2 )- The complete decay of a higher-order 
system to one binary and a number of single stars is expected 
to occur in less than 0.1 Myr. 

This raises the possibility that at least some classical and 
weak-lined T Tauri sources have been recently ejected from 
embedded multiple systems and are coeval with those sys- 
tems, not older than them. Ejected stars would leave a core in 
~ 100 kyr and depending on the size of a disc that they are 
able to retain during ejection would appear indistinguishable 
from classical and weak-lined T Tauri sources. Thus the classi- 
cal picture that protostellar sources evolve from class £>—>/—> 
// — » /// may need revision as it does not include (a) the possi- 
bility of creating later-type sources immediately through ejec- 
tions, and (b) the lack of a one-to-one relationship between the 
numbers of cores and the numbers of later-type sources. In par- 
ticular, the relative lifetime estimates of young stars assume a 



one-to-one correlation between embedded and non-embedded 
sources, the differences in relative numbers being due to the dif- 
ferent durations of each phase (e.g. Greene et al. 1994) which 
may be incorrect. 

2.2. Diluting the multiplicity fraction 

Ejections are predominantly of single stars (see Bate et al. 
2002, 2003; Delgado Donate et al. 2004; Goodwin et al. 
2004a,b for examples), and therefore ejections dilute / mu i t (e.g. 
Bouvier et al. 2001). A cluster which is (initially) composed 
of triple systems would have a / mmt of unity, however after the 
decay of these systems 1/3 of stars will be single, and 2/3 in bi- 
naries reducing f mu \ t to 0.5. Similarly, a population of quadru- 
ples will decay within less than 10 5 yr to a population with 
/mult = 1/3, assuming each quadruple decays to two single stars 
and one binary. 

The assumption is made here that the distances between 
members of the multiple system are approximately equal and 
that all of the members form at the same time. Simulations of 
multiple system formation within cores show that often hierar- 
chical systems can form with the members appearing at slightly 
different times (see Delgado Donate et al. 2004; Goodwin et 
al. 2004a,b). These multiples remain stable for at least 10 Myr 
(Delgado Donate et al. 2004). However, higher-order multiples 
in these simulations generally form as the result of the decay of 
even higher-order systems (N > 5 stars), so that / mmt is low for 
the whole population. 

In contrast to the above expectations from dynamical de- 
cay, star-forming regions are typically observed to have an ex- 
tremely high /muit (see Section 1). In addition to finding a sig- 
nificantly high /muit, Duchene et al. (2004) find no significant 
difference in / mu i t and the separation distribution between the 
physically different star-forming regions p Oph and Taurus, but 
a marked difference to the Galactic-field sample. Duchene et 
al. (2004; also see Haisch et al. 2004) find that f mu \ t is sig- 
nificantly higher in young embedded sources (38 + 8%) than 
in older flat-spectrum sources (22 + 3%) in both Taurus and 
p Oph (in both cases this is significantly higher than the field 
/muit of 14% for the same separation range). This suggests that 
there has been some dynamical decay of young higher-order 
multiple systems on a time-scale < 1 Myr. This observation 
also suggests that some of the late-type sources are single stars 
ejected from early-type multiple systems. 

In a cluster such as the Orion Nebula Cluster which has 
a /muit ~ 0.6 very similar to the field (Duchene 1999), about 
one star could have been ejected per binary system, i.e. cores 
may typically produce unstable triple systems, or alternatively, 
that roughly 1/3 cores only produced a single star while 2/3 
produced binaries. This leads to a basic constraint on star for- 
mation within cores in Orion of 1 < N < 3. On the other hand, 
without the need to invoke a different population of cores that 
form different numbers of objects to other star-forming regions, 
many (wide) primordial binaries can be destroyed by dynam- 
ical interactions in such a cluster (Kroupa et al. 2001). It is 
thus possible to match the Orion and the Galactic field / mmt 
and separation distribution if the initial population is comprised 
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mainly of binary systems with the T Tauri separation distribu- 
tion (Kroupa 1995a,b). 

It may be that ejected stars have escaped from the clus- 
ter (or at least travelled far enough not to be included in sam- 
ples) and so observations are of preferentially high f mu \ t re- 
gions. However, the ejection velocity of stars is ~ 1 km s - ' 
(Bate et al. 2002,2003; Delgado Donate et al. 2004; Goodwin 
et al. 2004a) which is comparable with the velocity dispersion 
of star-forming cores. Hence star-forming regions would be ex- 
pected to be mixed. Also, the potential well of larger clusters 
(such as Orion) is too deep to allow the loss of stars ejected at 
such velocities. We conclude that ejected stars and their pro- 
genitor systems should normally be coincident. 

The generally high / mu i t for pre-main sequence late-type 
stars uncovers an elementary discrepancy if cloud cores pro- 
duce N > 3 stars. 

2.3. Primordial triples 

A non-negligible fraction of T Tauri stars in the surveyed re- 
gions do appear to be hierarchical triples, probably > 25%. 
Leinert et al. (1993) find 4 higher-order multiples from 45 sys- 
tems in Taurus, at least 7 more are known with possibly more 
currently undetected. More recently, Koresko (2002) has used 
speckle holography at 2fim at Keck to search for tertiary com- 
panions among 14 pre-main sequence visual binary systems 
finding that possibly up to half of the imaged objects may be 
triples, although clear identification of tertiary companions was 
possible only in two of the possible 7 cases. Noticeably, a large 
fraction of Koresko's sample lies in the Ophiuchus cloud, a 
region that may be seen as a moderately rich cluster, clearly 
different from the Taurus cloud surveyed by Leinert et al., sug- 
gesting that this is a relatively general feature and not a "sin- 
gle cloud" phenomenon. In yet more dispersed star-forming 
regions (the TW Hya and MBM 12 associations), Brandeker 
et al. (2003) found a similarly high ratio of triples-to-binaries, 
though with limited statistical significance. 

In the Galactic field, Tokovinin & Smekhov (2002) suggest 
that 20-30% of multiple systems are higher-order hierarchical 
systems, consistent with the fraction of T Tauri stars suggesting 
that many of these young systems are stable as found in the 
simulations of Delgado Donate et al. (2004). The fraction of 
triple and quadruple systems per star, T — (T + Q)/(S + 2B + 
3T + 40 = 0.05 for the Galactic field, this corresponds to a 
higher-order multiplicity fraction of (T+Q)/(B+T + Q) = 18% 
of all multiple systems. 

The very high frequency of primordial binary and triple 
systems in the Taurus and p Oph star-forming regions suggest 
that very few ejections from multiple systems have occurred. 
This implies that stars must form in groups of 2 or 3 within 
cores. 

2.4. Hardening and close binaries 

Further support for our conclusion that 2 < N < 3 is provided 
by the separation distribution of systems. There are many bina- 
ries in Taurus with wide separations > 100 au (roughly double 



the number found in the field) which also implies that few ejec- 
tions can have occurred, as ejections rapidly harden binaries so 
that their semi-major axes are < 10 au - from an initial sepa- 
ration of 100 au, 2 or 3 ejections are sufficient (Anosova 1986; 
Umbreit et al. 2005). 

Ejections are expected to occur during the main (class 
0) accretion phase from basic decay-time arguments (see 
above), and also as seen in star formation simulations (Bate 
et al. 2002,2003; Delgado Donate et al. 2004; Goodwin et al. 
2004a,b; see also Reipurth & Clarke 2001). Ejections during 
this phase are far more effective at hardening binary systems 
(see Umbreit et al. 2005). This occurs as ejections harden a 
lower-mass system than the final system thus requiring less en- 
ergy. In addition, early hardening tends to lead to a popula- 
tion of equal-mass binaries, as the close low-mass component 
is more able to accrete high angular momentum material from 
the disc thus pushing the mass ratio towards unity (Whitworth 
et al. 1995; Bate & Bonnell 1997). Such systems are not ob- 
served in significant numbers in the field (Duquennoy & Mayor 
1991; Mazehetal. 1992). 

It is not expected that stars would form much closer to- 
gether than ~ 20 - 50 au (the separation between two minimum 
Jeans mass objects at ~ p cr \t, see Eqn. 1). Ejections are thought 
to be the primary method by which these systems are hardened. 
However, these systems are unlikely to have been hardened 
during the main accretion phase as they do not strongly favour 
equal-mass binaries, although short-period binaries do favour 
more similar component masses than wider binaries (Mazeh et 
al. 1992). They may however form from the (fairly late) decay 
of some fraction of the triple systems. 

Thus, in order to prevent a large population of close binaries 
being formed from the decay of unstable triple systems, N - 2 
must be a more common outcome of star formation than N = 3. 



2.5. Fragmentation consistent with data 

We propose the following picture that is consistent with the 
data. Of every 100 star forming cores within a cloud roughly 
60 will form binary and 40 will form triple systems. Of those 
which form triples, around 25 are long-lived hierarchical sys- 
tems, whilst the other 15 triples decay to 15 binaries and 15 
single low-mass stars and/or brown dwarfs. The 15 binary sys- 
tems that result from dynamical decay form the close binary 
population. 

This leads to 115 systems, 15 of which are single stars, 75 
are binaries and 25 triples with / mu i t = 100/1 15 = 0.87 for the 
whole population. As the ejected stars were the lowest mass 
members of the triple systems, / mu i t remains ~ 1 for larger 
> 0.3-0.5M© stars, but is low among low-mass stars and brown 
dwarfs. Such a primordial population would then either remain 
essentially unevolved in Taurus-like star-forming regions, or 
would evolve dynamically in dense Orion-like environments 
on a cluster crossing-time scale (Kroupa 2000; Kroupa et al. 
2001; Kroupa et al. 2003). 
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2.6. Brown dwarfs 

The dynamical ejection of sub-stellar embryos has been sug- 
gested as a formation mechanism for brown dwarfs (Reipurth 
2000; Reipurth & Clarke 2001). Ejections are expected to oc- 
cur during the main accretion phase (see above), indeed brown 
dwarfs must be ejected shortly after they form as the accretion 
rates of very young sources rapidly grow the embryos to stel- 
lar masses (which has been suggested as an explanation of the 
brown dwarf desert, Reipurth & Clarke 2003). The ejection of 
brown dwarfs, despite their low mass, would still have a signif- 
icant effect on the separation (and hence mass ratios) of young 
multiple systems if they occur when all of the components have 
not accreted a significant amount of mass (see also Umbreit et 
al. 2005). Thus, ejections of brown dwarfs during the main ac- 
cretion and star formation phase 1 must be limited in order to 
avoid creating a large population of close, equal-mass binaries. 

The available star-count data in the Taurus star-forming 
region and in Orion constrains this process to about one 
brown dwarf being ejected per four stars in both environments 
(Kroupa & Bouvier 2003). So this scenario is not ruled-out 
as a brown dwarf formation mechanism, especially if embryo 
formation can be delayed and brown dwarfs ejected from sys- 
tems whose other members have a significantly higher mass, 
in which case the hardening of the parent system would be re- 
duced. It should be emphasised that, even if all brown dwarfs 
are ejected embryos, the limit of one brown dwarf per 4 stars 
(or per two stellar binary systems) still constrains the total num- 
ber of ejections to be low. 

2. 7. Cluster remnants as a source of high-order 
systems 

As stated above, many of the observed N = 3 systems in the 
field may be primordial, i.e. may have formed as stable hier- 
archical systems from molecular cloud cores. However, some 
triple and perhaps most higher-order multiple systems may also 
be formed by star cluster remnants. 

de La Fuente Marcos (1997, 1998) points out that cluster 
remnants are characterised by a few binaries, ie. highly hierar- 
chical high-order multiples, and that the number of such objects 
can be large given that the majority of Galactic-field stars form 
in modest clusters. Evaporating clusters decay into highly hier- 
archical systems of N < 15 members. These remnants will pre- 
sumably decay further over time possibly producing quadruple 
or quintuple systems. We can estimate the number of highly 
hierarchical N > 3 systems under this scenario as follows: 

For a power-law embedded cluster mass function, £ e ci = 
k M^ v where £ e d dM ec i is the number of embedded clusters 
with mass in the interval M ec i,M ec i + dM ec i and M ec i is the 
mass in embedded cluster stars, the number of hierarchical 
high-order stellar systems is equal to the number of embed- 
ded clusters in a mass interval M ec i >m j n to M ec i jmax , N ec i = 
since each cluster leaves one such system. 

1 It is expected that star formation would become difficult after the 
class phase as most of the gas available for star formation has already 
been accreted. 



The number of stars formed in the embedded cluster population 
during some epoch is N st = J M ^\M ed /m m ) £ e d dM ed , where 
m av = 0.4 M is the average stellar mass which we take to be 
constant for the present estimate. Thus, for M ecl min = 5 M Q (a 
cluster of a dozen stars, Taurus-like object, the least-massive 
star-forming "unit" known) and M ec i imax = 500 M (since most 
Galactic-field stars appear to stem from rather modest clusters, 
Kroupa 1995a; Adams & Myers 2001) and for /? = 2 (Lada 
& Lada 2003) we have N eci = 0.19 and N st = 4.43 so that the 
above defined quantity f = N ec \/N st ~ 4%. Therefore the num- 
ber of cluster remnants may be completely sufficient to account 
for all observed higher-order multiple systems. 

3. Conclusions 

The high multiplicity fraction / mu i t of pre-main sequence stars 
implies that most stars form in multiple systems within dense 
molecular cores. The decay of young multiple systems gener- 
ally within < 0.1 Myr has two main effects: diluting the mul- 
tiple population through the ejection of predominantly single 
stars, and hardening the remaining systems. 

The observed large fraction of ~ 10 6 yr old, wide binary 
and hierarchical triple systems implies that (a) not much dilu- 
tion of the binary population can have occurred through ejec- 
tions, and (b) few ejections can have occurred as systems have 
not been significantly hardened. 

This suggests that most cores form only 2 < N < 3 stars. 
Typical numbers that are consistent with observational con- 
straints are that of 100 birth systems 40 are triple and 60 bi- 
naries. Of the 40 triples 25 are long-lived hierarchies while 15 
decay to 15 close binaries and 15 single stars. This is then the 
primordial population as seen dynamically unevolved in Taurus 
and p Oph (Kroupa et al. 2003). The destruction of binary sys- 
tems by dynamical interactions in dense clusters such as the 
Orion Nebula cluster further dilutes the binary fraction (Kroupa 
et al. 2001). Wide and low-mass binaries are preferentially de- 
stroyed, and the majority of single low-mass stars can be ex- 
plained through this process and do not require ejections. Such 
dynamical disruptions occur in the early stages of cluster evo- 
lution and explain the rapid drop in the binary fraction observed 
between T Tauri stars and older clusters and the Galactic-field 
main sequence population. The observed variations may be en- 
tirely understood through dynamical disruption of wide binary 
systems in differently crowded regions. 

This would imply a remarkable invariance of core fragmen- 
tation to differing star-formation conditions. Motte & Andre 
(2001) have shown that protostellar cores show different sizes, 
densities and edge sharpness if they are located in cluster- 
like environments or sparser associations. Nevertheless, the re- 
markable success in describing a variety of stellar populations 
with the same primordial binary population does suggest that 
core fragmentation may indeed be about invariant on the scale 
of a few hundred au or less. Such an invariance may not be 
unexpected as, whilst large-scale core properties vary signifi- 
cantly, it is only after collapse through several orders of mag- 
nitude in density that stars form, and such a collapse may make 
many large-scale properties irrelevant. 
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Ejections from cloud cores can lead to brown-dwarfs as 
ejected embryos, but the frequency of such events is already 
limited by available data to not more than 1 brown dwarf per 
two stellar binaries (or one per four stars). 

Cluster remnants can form highly hierarchical high-order 
multiple systems. These contribute to the observed number of 
N > 3 stellar systems. Cluster remnants may account for all 
observed higher-order (N > 3) multiple systems. Therefore, 
possibly only a small fraction of molecular cloud cores have 
to be able to form long-lived stable hierarchical higher-order 
multiple systems. 

These conclusions place strong constraints on theories of 
star formation. For any theory of star formation to match ob- 
servations the majority of cores must fragment into multiple 
objects. However, they can usually only fragment into 2 or 3 
stars. The currently available theoretical results appear to be 
inconsistent with this, as the cloud-core fragmentation calcula- 
tions typically form N = 5 - 10 fragments per core. This dis- 
crepancy may not be unexpected because the theoretical work 
misses the very essential physics of stellar feedback, and can- 
not take into account magnetic forces that may play a role in the 
fragmentation of a core, and simplifies the complex and impor- 
tant thermal physics (in particular simplifying the equation of 
state of the gas). 
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